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ABSTRACT 



Context. Large amplitude oscillations of solar filaments is a phenomenon known for more than half a century Recently, 
a new mode of oscillations, characterized by periodical plasma motions along the filament axis, was discovered. 
Aims. We analyze such an event, recorded on 23 January 2002 in Big Bear Solar Observatory Ha filtergrams, in order 
to infer the triggering mechanism and the nature of the restoring force. 

Methods. Motion along the filament axis of a distinct buldge-like feature was traced, to quantify the kinematics of the 
oscillatory motion. The data were fitted by a damped sine function, to estimate the basic parameters of the oscillations. 
In order to identify the triggering mechanism, morphological changes in the vicinity of the filament were analyzed. 
Results. The observed oscillations of the plasma along the filament was characterized by an initial displacement of 24 Mm, 
initial velocity amplitude of 51 kms~^, period of 50 min, and damping time of 115 min. We interpret the trigger in 
terms of poloidal magnetic flux injection by magnetic reconnection at one of the filament legs. The restoring force is 
caused by the magnetic pressure gradient along the filament axis. The period of oscillations, derived from the linearized 
equation of motion (harmonic oscillator) can be expressed as P = 7tV2L/va^ ~ 4AL/va^j where va^ = B^o/yf/Iop 
represents the Alfven speed based on the equilibrium poloidal field B^o. 

Conclusions. Combination of our measurements with some previous observations of the same kind of oscillations shows 
a good agreement with the proposed interpretation. 
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1. Introduction 

Solar coronal structures are frequently subject to oscil- 
lations of various modes and time/spatial scales (e.g., 
Aschwanden I2003p . The longest known phenomenon of 
this kind are so called winking filaments (cf., Ramsey & 
Smith . 1966 ). They represent large-amplitude large-scale os- 
. ^ , dilations of prominences observed on the solar disc, most 
often triggered by disturbances coming from distant flares. 
^ Later on, various modalities of prominence oscillatory mo- 
^ tions were reported (e.g., Kleczek & Kuperus 1969, Malville 
" " " & Schindler HnHIl Vrsnak[198l Wiehret al. .l984 Ts ubaki 
& Takenchi [198 1 Vrs nak et al. [T99Q1 Jing et al. l2QQ3l 
Isobe & Tripathi 2006) . The phenomenon is generally in- 
terpreted in terms of different magnetohydrodynamical 
(MHD) wave modes (for a classification we refer to re- 
views by Tsubaki |1988t Vrsnak d993| Roberts i2000| Oliver 
& Ballester 2002]). 

Most of reported large-amplitude large-scale oscillations 
happen perpendicular to the prominence axis. The restor- 
ing force in this type of oscillations generally could be ex- 
plained in terms of the magnetic tension (e.g., Hvder [T966| 
Kleczek & Kuperus 1969, Vrsnak 1984, Vrsnak et al. [15501) . 
On the other hand, the oscillation damping in the corona 
is attributed either to the "aerodynamic" drag (i.e., the 
energy loss by emission of waves into the ambient corona; 
e.g., Kleczek & Kuperus 1969), or to various dissipative pro- 
cesses (e.g., Hyder 1966, Nakariakov et al. 1999, Terradas 
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et al. [200n Ofman & Aschwanden \2002\ Verwichte et 

Recently, Jing et al. (120031 12006|) reported a new mode 
of oscillations in filaments/prominences, where the oscilla- 
tory motion happens along the prominence axis. In this pa- 
per we report another example observed in Ha filtergrams 
and propose an explanation for the triggering process and 
the restoring force. 



2. Observations 

The filament that we discuss was located near the center 
of the solar disk on 2002 January 23 and was observed in 
full-disk Ha filtergrams at I'Ypixel resolution with a time 
cadence of 1 min at the Big Bear Solar Observatory (BBSO) 
Global Ha Network Station (Steinegger et al. 2000). The 
accompanying movie shows the evolution of the oscillating 
filament between 17 and 22 UT; snapshots are presented in 
Fig. [TJ During this time span, the filament was activated 
and showed oscillations along its main axis which could be 
followed for 5 cycles. 

Before the measurements, we co-aligned all the images 
to the reference image taken at 17:00 UT (hereinafter used 
as t = 0), correcting for solar differential rotation and also 
for image jittering by the means of cross-correlation tech- 
niques. For measuring the oscillatory motions of the fila- 
ment, we did not use all the available images but on aver- 
age each third image of the series, giving an average time 
cadence of 3 min. 



2 



B. Vrsnak et al.: Large amplitude oscillatory motion along a solar filament 



The main part of the filament, 110-120 Mm long, was 
following the straight magnetic inversion line oriented in 
Southeast-Northwest direction. Its Northwest leg, a curved 
thread of length 20 - 30 Mm, was rooted in a small active 
region, NOAA 9793 (Fig. [2]). The Southeast leg was rooted 
in the quiet chromosphere. For the total length of the fila- 
ment we take 2L = 140 Mm. 

The kinematics of the observed oscillations was deter- 
mined by visually inspecting the Ra image sequence, and 
measuring the displacement along the filament axis of the 
most prominent (dark) bulge-like feature in the filament 
(Fig. [U see also the accompanying movie). The measured 
displacements are presented in Fig. [3^, and the correspond- 
ing velocities in Fig. [3)3. Note that plasma motions could 
be seen all along the filament, with the amplitude decreas- 
ing towards the filament legs, i.e., the oscillation amplitude 
was largest in the middle part of the filament, where the 
bulge-like feature was located. A mechanical analogue of 
the observed oscillation would be, e.g., a longitudinal- mode 
standing wave on a slinky spring fixed at both ends. 

It seems that the filament oscillation was a response to 
an energy release in the active region, where a weak flare- 
like brightening appeared in Ha at 17:33 UT in the region 
of the Northwest footpoint of the filament (see the second 
panel in the first row of Fig. (TJ see also the accompany- 
ing movie). The Ha sub flare achieved a maximum around 
17:45 UT, as can be seen in Fig. [Sb where we show the 
Ha light curve, measured at the western kernel of the sub- 
flare (see the inset in Fig. [3t; the eastern kernel was partly 
obscured by the activated filament, thus it was not possi- 
ble to measure the light curve in this part of the subflare. 
The energy release associated with the Ha brightening was 
obviously weak, since no increase in the GOES full-disk 
integrated soft X-ray flux was observed. 

The subflare caused the filament activation, where var- 
ious forms of motions and restructuring could be observed. 
Due to complex motions and morphological changes of the 
filament, in the period 17:35-18:10 it was practically im- 
possible to identify, and measure reliably, the position of 
the filament feature which was later on used to trace the 
oscillations. In the period between 18:00 and 18:10 UT the 
motions gradually became more ordered, mainly showing a 
flow of the filament plasma in the Southeast direction, away 
from the active region. After this initial plasma displace- 
ment, the first Northwest directed swing started around 
18:10 UT (Fig. [3^; the first oscillation measurement was 
performed at 08:12 UT). 

The measured displacement curve, d{t)^ of the oscilla- 
tory filament motion was fitted by a damped sine function 
of the form 



d{t) = xq sin 



27rt 



(f I exp 



do 



(1) 



with amplitude xq, period P, phase (/p, damping time r 
and equilibrium position do. Fig. [3^ shows the observed 
oscillatory motion of the filament (with do already sub- 
tracted), together with the fitted curve. In Fig. [3)3 we plot 
the corresponding velocities obtained by numerical differen- 
tiation of the displacement data points (applying a 3-point, 
Lagrangian interpolation) together with the time derivative 
of the d{t) fit. 

The basic parameters of the oscillations are: 

— period P = 49.8 min; 
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Fig. 3. a) Displacement measurements of the oscillating fil- 
ament together with the applied fit. b) Velocities derived 
via numerical differentiation from the displacement mea- 
surements (for the procedure see main text) shown together 
with the time- derivative of the d(i) fit shown in panel a), 
c) Ha light curve measured in the western footpoint of the 
associated subflare indicated by the box in the inserted im- 
age. 



initial amplitude xo = 24.1 Mm; 

initial velocity amplitude i;o = 51 kms~^; 

damping time r = 114.9 min. 



3. Interpretation and discussion 

The oscillations described by Eq. (1) are the solution of the 
equation of motion of the form: 



- cJqX = . 



(2) 



The cyclic frequency of oscillations a; is related to the fre- 
quency of free oscillations ujo and the damping rate b by 
the expression uP' = uj^ — . From the period of oscilla- 
tions and the damping time determined in Sect. 2 one finds 
the cyclic frequency of oscillations and the damping rate as 
uj = 2tt/P = 2.1 X 10-^ s-i and S = l/r = 1.5 x IQ-"^ s-\ 
respectively. Since co ^ it can be concluded that cjq ^ ^, 
so in the following we take for the eigenmode frequency of 
the system the value cjq ^ 2 x 10~^ s~^. From the cyclic fre- 
quency ooo and the amplitude xo one finds the initial veloc- 
ity and acceleration amplitudes, vo — oooxo = 50.8 kms~^ 
and ao = oOqXo = 107 ms~^. 
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Fig. 1. Sequence of BBSO Ha images for a 210^^ x 210^^ subfield. For each image, the observing time (in minutes after 
17:00 UT) is annotated. The arrows connect the determined displacement for consecutive images taken around the 
oscillation peaks for each of the five cycles (C1-C5). The arrowhead indicates the determined center of the followed 
filament feature in the actual image, the circle drawn at the other end of the arrow indicates the location of the same 
feature as determined from the previous image shown in the figure. Thus, the decreasing lengths of the arrows directly 
reflect the instantaneous oscillation amplitudes and their damping. See also the accompanying movie. 
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Fig. 2. MDI/SOHO continuum image, MDI/SOHO magnetogram and BBSO Ha filtergram at 17 UT of active regions 
NOAA 9787 and NOAA 9793 with the filament. Units of the x- and y-axis are in arc sec. 



The estimated value of acceleration gives a constraint 
to the nature of the restoring force. Let us first assume that 
the restoring force is governed by the gas pressure gradi- 
ent. It can be approximately estimated as Ap/2L, where 
we assume that there is a pressure difference over the 
filament length 2L. Such a pressure gradient would give rise 
to an acceleration of a = Ap/2Lp, where p is the promi- 
nence density. Writing down an order of magnitude esti- 



mate 2La = Ap/p ^ p/p ^ c^, one finds that the sound 
speed has to be Cg = \plLa > 100 kms~^. Such a sound 
speed requires temperatures > 10^ K, which is typical for 
the coronal plasma. So, the required sound speed is an or- 
der of magnitude too high, since the temperature of the 
prominence plasma is ^ 10"^ K, corresponding to a sound 
speed in the order of Cs ^ 10 kms~^. A similar conclu- 
sion could be drawn by taking into account that the pe- 
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riod of oscillation can be roughly identified with the sound 
wave travel time along the filament axis, tg ^ 2L/cs^ i.e., 
expressing the period as P 2ts = ^L/cs. Substituting 
2L = 140 Mm, we find again Cg > 100 kms~^ (note that 
using the oscillation amplitude xq instead of the filament 
length 21/ would give an estimate of the velocity amplitude 
and not the oscillation-related wave velocity, i.e., the sound 
speed). 

So, the pressure pulse scenario would be possible only if 
the pressure of the filament plasma was increased by a fac- 
tor of 100. That corresponds to a temperature increase from 
To ~ 10^ K to T > 10^ K, but no signature of such plasma 
heating was observed in the TRACE 171 A (T - 1 MK) 
EUV images. Another possibility is that the restoring force 
is of magnetic origin. For example, let us consider that the 
filament is embedded in a flux rope and that at a certain 
moment an additional poloidal flux is injected into the rope 
at one of its legs. The flux could be injected from below the 
photosphere, or by magnetic field reconnection (see, e.g., 
Uchida et al. "2001, Jibben & Canfield 2004). Indeed, Ha 
fiare-like brightenings appeared just before/around the on- 
set of oscillations in the region of the Northwest footpoint 
of the filament (see Sect. 2), which is likely to be a signature 
of reconnection. On the other hand, no sign of the emerg- 
ing fiux process could be detected in the region around 
the filament in the SOHO/MDI magnetograms which were 
available at a time cadence of 1 min. 

The excess poloidal field creates a magnetic pressure 
gradient and a torque, driving a combined axial and ro- 
tational motion of plasma (e.g., Jockers Il978t Uchida et 
al. 120011) . If the injection lasts shorter than the Alfven travel 
time along the flux rope axis, it causes a perturbation that 
propagates along the rope (Uchida et al. [200 ID and evolves 



X = x/L, we get the equation of motion of the harmonic 
oscillator: 



into an oscillatory motion along its axis. Here we assume 
that the wave caused by the injection is reflected at the op- 
posite footpoint of the rope, and that the damping rate is 
smaller than the oscillation cyclic frequency. Also note that 
since Cg < va the plasma is carried by B^p since plasma 
motions along the field lines are slower than translation of 
poloidal field (fast-mode MHD wave). 

Let us consider a flux rope of a constant width, where 
a certain amount of poloidal flux is injected at one of its 
footpoints (similar to that shown in Fig. 5 of Jibben & 
Canfield 2004), over an interval shorter than the eigenmode 
period of the related oscillation mode. We divide the flux 
rope in equilibrium into two segments of length L, symmet- 
ric with respect to the flux rope midpoint. In the simplest 
form, the average poloidal magnetic field in the two seg- 
ments, after perturbing the flux rope, can be expressed as 
B^i = B^qL/{L - x) and B^2 = B^qL/{L^x). Here, B^q 
is the poloidal field in the equilibrium state and x is the 
longitudinal displacement of the plasma element separat- 
ing the two segments (in the equilibrium taken to be at 
the filament barycenter, located at the midpoint of the flux 
rope). We use these relations to express the corresponding 
magnetic pressures, Prn = ^^/^/io, which we substitute into 
the simplified equation of motion: 



px 



Apr, 

L 



^^0 

2/io 



1 



(L - x)2 (L + xy 



(3) 



where Apm/L represents approximately the gradient of the 
poloidal- field magnetic pressure. After linearizing Eq. (3), 
and expressing the displacement in the dimensionless form 



X 



L2 



X. 



(4) 



where va^p = B^q / yfjUop represents the Alfven speed based 
on the equilibrium poloidal field B^^q. 

Equation (4) implies loq = 2v\^/L'^^ and consequently, 
P = i^\f2LlvA^ ^ ^ALjvA^' Utilizing the observed pa- 
rameters, we find va^ ^ 100 kms~^. Assuming that the 
number density of the prominence plasma is in the range 
n = 10^^-10^^ cm~^ (cf. Tandberg-Hanssen 11995j), one 
finds 5(^0 ^5-15 gauss. 

In this respect we note the filament was composed of 
helical-like fine structure patterns, twisted around the fil- 
ament axis. Bearing in mind that the frozen-in condition 
is satisfied in the prominence plasma, such patterns are 
indicative of the flux-rope helical magnetic field (Vrsnak 
et al. IT99T]) . In the period 18:25-18:40 UT the internal 
structure was seen clearly enough to provide measurements 
of the pitch angle of helical threads. The measured val- 
ues of the pitch angle range between 20° and 30°, with 
the mean value tan'i^ = X = 24° ± 4°. Bearing in mind 
X — Bip/B\\^ where B^p and B\\ are the poloidal and axial 
magnetic field components (Vrsnak et al. ll991|) , and consid- 
ering the poloidal field of 5-15 gauss, one finds B\\ ^10- 
30 gauss, which is a reasonable value for quiescent promi- 
nences (Tandberg-Hanssen "1995) . 

According to Vrsnak (1990), who treated the stability of 
a semitoroidal flux rope anchored at both legs in the pho- 
tosphere, the prominence with pitch angle ^ = 20° - 30° 
should be stable, which is consistent with the fact that the 
filament was severely activated, but did not erupt. A similar 
conclusion can be drawn if the total number of turns of the 
helical field line is considered. Taking into account that the 
filament width was around 2r ~ 2 arcsec and employing the 
expression for the pitch- length A = 27rr/X, we find the to- 
tal number of turns of the helical field line TV = 2L/A ~ 1.4. 
This is again consistent with the model stability-criteria by 
Vrsnak [T990I (for the unstable cases of highly twisted quies- 
cent prominences see, e.g., Vrsnak et al. . 1991. or Karlicky & 
Simberova l2002p . However, it should be noted that some- 
times active region prominences erupt at smaller number 
of turns {N > 1; see, e.g.. Rust & Kumar 1996^ Regnier & 
Amari[2nQl Salman et al. [2QQ7|) . 

Regarding the P oc L relationship, we have to discuss 
Fig. 5a of Jing et al. (|2006p , where in the figure caption it is 
stated that the x-axis of the graph represents the "length" . 
However, checking Table I therein, one finds that in fact 
the X-axis represents the oscillation amplitude, which we 
previously denoted as xq. Consequently, the lack of corre- 
lation does not really mean that there is no P cx L rela- 
tionship. In this respect it is instructive to comment also 
Fig. 5b of Jing et al. (2006) which shows proportionality of 
the velocity amplitude vq and "length", the latter in fact 
representing the amplitude xq. Such a proportionality is ex- 
pected for any harmonic oscillator since vq = cooXq. Bearing 
in mind Eq. (4), i.e., uoq oc 1/L, one would expect an in- 
verse proportionality between vq and L in the case of equal 
amplitudes. 

Out of four events studied by Jing et al. ([20031 [2QQ6|) . 
in two cases it is possible to estimate the overall length 
of the oscillating filaments. From filtergrams displayed in 
Fig. 1 of Jing et al. (f2QQ3l) and Fig. 3 of Jing et al. (;20Q6|) . 
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showing filaments characterized by periods P = 80 and 
100 min, we estimated the filament lengths to approxi- 
mately 2L ^ 200 Mm and 2L ^ 300 Mm, respectively. 
Combining these values with our measurements (P = 
50 min; 2L ^ 140 Mm), we find that the data follow approx- 
imately the P (X L dependence. The linear least squares fit 
gives a slope of 0.6 min/Mm, or 0.7 min/Mm if the fit is 
fixed at the origin. If we compare these values with the 
previously derived relationship P ^ AAL/va^^ we find that 
these slopes correspond to VAif ~ 100-120 kms~^, which 
is in agreement with the value of VAif that we estimated 
earlier for our oscillating filament. 

Finally, let us note that the poloidal flux injection 
changes the poloidal-to- axial field ratio X = B^/B\\, which 
should cause also oscillations perpendicular to the promi- 
nence axis, since the equilibrium height of the prominence 
is directly related to the ratio X (e.g., Vrsnak [T984| IT990I) . 
Indeed, careful inspection of the filament motion reveals 
such transversal oscillations of a small amplitude, syn- 
chronous with the longitudinal ones (see the accompanying 
movie). The analysis of the interplay between the longi- 
tudinal, transversal, torsional (pitch- angle), and radial os- 
cillations requires a meticulous treatment, which will be 
presented in a separate paper. 



4. Conclusion 

The observed oscillations of the plasma along the filament 
was characterized by a period of 50 min, velocity ampli- 
tude of 50 kms~^, and damping time of 115 min. Our anal- 
ysis indicates that the oscillations of the filament plasma 
along its axis were driven magnetically as in other types of 
prominence oscillations. We propose that the oscillations 
were triggered by an injection of poloidal magnetic field 
into the flux rope, most probably by the reconnect ion as- 
sociated with a subflare that took place at the Northwest 
leg of the filament. The flux injection is expected to cause 
a translatory and torsional motion of plasma, propagating 
toward the other leg of the filament. After the reflection of 
the wave, the oscillations along the filament axis develop, 
showing an order of magnitude weaker damping than in 
perpendicular oscillations of filaments (for the latter see, 
e.g., Hyder 1966). Most probably, perpendicular oscillations 
decay much faster because they affect the ambient corona 
more than motions along the flux tube, i.e., the energy flux 
carried away by MHD waves is much larger. 

The period in our event is shorter than in analogous 
oscillations analyzed by Jing et al. (20031 [2QQ6|) , where the 
periods were in the range from 80 to 160 min. In two of 
these events the filament length could be measured, and 
when supplemented with our event, we find that the periods 
follow the P (X L relationship. The slope of the relationship 
is consistent with a poloidal magnetic field in the order of 
10 gauss. 

The damping time r is also shorter than in the events 
analyzed by Jing et al. (|2003[ l2QQ6|) . If our measurements 
would be inserted into Fig. 5c of Jing et al. (2006), the 
corresponding data-point would be located very close the 
lower end of the regression line presented therein. The linear 
least squares fit, with the axis intercept fixed at the origin, 
would read r = 3.5P. 
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Fig. 4. Solutions of Eq. (5) shown for amplitudes Xq = 
0.25, 0.5, 0.75, and 0.9 (dashed, dotted, dot-dashed and dot- 
dot-dashed line respectively), together with the harmonic 
oscillator solution for Xq = 0.25 (full line). 

Appendix: The nonlinearity aspect 

Since the filament oscillations were characterized by large 
amplitudes, it is instructive to check how much the 
harmonic-oscillator approximation differs from the solu- 
tions of Eq. (3). After introducing the dimensionless dis- 
placement X = x/L^ Eq. (3) can be rewritten in the form: 



where uJq = 2v\^/L'^. In Fig. [4] we show the solutions of 
Eq. (5) for amplitudes Xq = 0.25, 0.5, 0.75, and 0.9 to- 
gether with the harmonic oscillator solution for Xq = 0.25. 
We have chosen cjq = 4 x 10~^ s~^, approximately cor- 
responding to P = 50 min. Inspecting Fig. [H one finds 
that the period decreases with increasing amplitude: it is 
about 2 % shorter than the harmonic oscillator period for 
Xo = 0.25, around 10% shorter at Xq = 0.5, and 25% 
at Xq = 0.75. The deviation from the harmonic oscillator 
period at Xq = 0.9 becomes larger than 40%. Bearing in 
mind the damping, the increase of the period with decreas- 
ing amplitude means that the period should increase as the 
oscillations attenuate. Indeed, going back to Fig. O we find 
that the measurements might indicate such a trend towards 
the end of oscillations, in particular after t > 250 min. 
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